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Abstract. Rigel is a bright nearby B8 la supergiant star. We observed Rigel for 27.7 
days with the MOST satellite and monitored it's optical spectrum for over 6 years. Radial 
velocity (RV hereafter) curve manifestly calls for tens of modes for prewhitening. We 
conservatively report on the first 19 significant modes to avoid misdetection of aliases. The 
variability periods range from about 75 days down to a day. All radial modes are stable. 
We speculate the presence of gravity-dominated mixed-modes excited by e— mechanism from 
the Hydrogen burning shell on top of He burning core. 



1. Introduction 

Blue supergiant phase is a very transient, short-lived evolutionary phase of massive stars (M > 
15M Q ), and hence the study of the internal structure of these core-collapse supernova progenitors 
could bring important constraints on the mass, size, composition and rotation rate of the He 
core prior to the collapse. This, however, depends critically on how different mixing mechanisms 
such as overshooting, semi-convection and rotation are treated during the main sequence (MS) 
phase of massive stars (Heger et al. 2000; Noels et al. 2010). In this respect, we find Rigel an 
excellent testbed for such investigations. 

Saio et al. (2006) detected 48 pulsation modes in the MOST post-MS target HD 163899 (B2 
Ib/II). They explained the excitation of observed g- modes in terms of mode reflection from the 
deep lying intermediate convective zone (ICZ) when they arrive at the appropriate phase to the 
edge of the ICZ. Godart et al. (2009) theoretically investigated the destructive effects of mass 
loss and core overshooting on the location and extent of the ICZ, and the asteroseismic potential 
of this class of stars. Since a limited number of modes can survive the huge radiative damping 
in the core, they propose to artificially impose the inner boundary condition at the base of ICZ. 
According to Saio (2011) expected instabilities for stars like Rigel are radial modes in addition 
to / = 1 and 2 convection g~ modes. However, Gautschy (2009) predicts no instability for stars 
with 3.95 < logoff < 4.15 (Rigel lies in this gap). This is why asteroseismology of Rigel could 
be model-dependent and demanding. 

Several observed quantities of Rigel which impose valuable constraints on the equilibrium 
model of the star are already measured. They are the Hipparcos parallax ir = 3.79 ± 0.34 
mas, effective temperature T g = 12100 ± 150 K, surface gravity log g = 1.75 ± 0.10, luminosity 
log(L/L©) ~ 5.07 ± 0.10, metalicity [M/H]=-0.06 ± 0.10 (Przybilla et al. 2006), surface He 
abundance Y s = 0.32 ± 0.04, vsini ~ 25 ± 3 km s _1 (Przybilla et al. 2010; Simon-Diaz et al. 
2010), and limb darkened angular diameter #ld = 2.75 ±0.01 mas. Moreover, Rigel's variability 
in the Ha line profile (Kaufer et al. 1996; Morisson et al. 2008) and other metal-line equivalent 
widths due to non-radial pulsations (Kaufer et al. 1997) and mass loss (M = 1.5 ±0.4 x 10~ M & 
yr , Chesneau et al. 2010) are already published. These observed values are later tried to be 
reproduced when modeling the star. When modeling the equilibrium state of Rigel (section 3.) 
we attempted to match the above mentioned observed quantities with the modeled quantities. 



1 



2 



E. Moravveji et al. 



2. Variability Time Scales Prom Radial Velocity Monitoring 

Photometry of Rigel was carried out with MOST (Walker et al. 2003) satellite for 27.7 days 
from 15 November to 13 December 2009. Since our photometry time-baseline is of the same 
order of the long-period modes no definite frequencies are found in Lomb-Scargle periodogram. 
The RV monitoring, in return, extends over 6 years but suffers from daily and annual gaps. We 
detected 19 significant harmonics using SigSpec software (Reegen 2007) above sig=15 threshold; 
the results are in complete agreement with those of Perod04 software (Lenz k, Breger 2005) with 
SNR>4.6. The extensive list of harmonics will be presented in Moravveji et al. (in preparation). 
The two lowest frequencies occur at 0.0152 and 0.0134 d" 1 with semi-amplitudes of 1.126 and 
0.838 km s -1 , respectively, and the highest frequency occurs at 0.8203 d _1 . 

3. Equilibrium Modeling and Non-adiabatic Radial and Non-radial Pulsations 

Stellar structure and evolution modeling which includes rotational mixing (Heger et al. 2000), 
mass loss, and element diffusion subject to constraints listed in section 1. is carried out with the 
MESA code (Paxton et al. 2011). Schwarzschild criterion for convection with «mlt = 1-6 is used. 
So, the main source of energy is a slight He burning in the core followed by H burning in shell. 
But models with Ledoux convective criterion have no He burning core (Godart et al. 2009). 

Since non-adiabaticity can play a critical role in excitation and damping of modes in massive 
stars, we treat modal stability analysis employing Granada Oscillatin Code (GraCo, Moya et al. 
2004). For the input equilibrium model from MESA, first we scan for < I < 3 modes in 
adiabatic frame, and next we solve the full non-adiabatic equations. Even without imposing 
the inner boundary condition at the base of the ICZ (suggested by Godart et al. 2009), there 
are several unstable high order low-degree mixed gravity modes as eigensolutions. Fundamental 
radial mode and its overtones are all stable. The radiative layer below the ICZ where H- 
burning takes place shows significant contribution in eigenfunction of excited modes. After 
analyzing the work integral W and its derivative dW/dhiT, we speculate that the pulsation 
driving arises from the e— mechanism (Unno et al. 1989). This is powered by the sensitivity of 
the CNO burning network to the temperature perturbations below the ICZ of logT ~ 7.7 which 
dominates over k— mechanism around Iron bump of opacity located at logT pa 5.2. If securely 
proven, this approach will open a venue towards studies of near-core condition in progenitors of 
SN II progenitors. This part of the study is still under progress. 
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